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I.  INTRODUCTION 


A.  Background 

It  has  been  known  for  some  time  that  spin-stabilized  projectiles  with 
liquid  payloads  can  suffer  flight  instabilities  related  to  the  presence  of  the 
liquid1.  It  has  also  been  known  that  a  rotating  fluid  can  sustain  inertial 
oscillations2,  and  that  these  oscillations  are  constrained  to  certain  discrete 
natural  frequencies  if  the  fluid  is  contained.  In  1959,  Stewartson3 
demonstrated  that  a  rotating  rigid  body  having  a  right  circular  cylindrical 
cavity  filled  with  liquid  can  suffer  an  instability  of  its  fast  precessional 
motion,  if  the  frequency  of  the  body  is  sufficiently  close  to  one  of  the 
natural  frequencies  of  a  nonsymmetric  mode  of  inertial  oscillation  of  the 
liquid.  Stewartson  and,  independently,  Phillips4,  calculated  the  eigen- 
frequencies  for  an  inviscid  liquid,  and  Stewartson  presented  a  linear  theory 
giving  the  frequency  of  precession  of  the  rotating  body-liquid  system,  and  the 
rate  at  which  the  angle  of  precession  grows  or  decays  with  time.  Wedemeyer5’6 
used  boundary  layer-type  corrections  to  determine  the  change  in  eigenfre- 
quencies  from  inviscid  values  caused  by  small  viscosity,  and  proposed  that  the 
change  in  the  frequency  and  growth  rate  of  system  precession  due  to  small 
viscosity  could  be  entirely  explained  (to  first  order  accuracy  in  the  Ekman 
number)  by  the  change  in  liquid  eigenfrequencies.  Stewartson' s  theory  with 
Wedemeyer' s  extension,  hereafter  referred  to  as  the  Stewartson-Wedemeyer 
theory,  comprises  the  theoretical  basis  for  the  principal  tool1  now  available 
for  the  design  of  stable  liquid-filled  projectiles.  Important  assumptions 
made  in  the  derivation  of  the  theory  which  limit  its  range  of  applicability 
are  that  (i)  the  angle  of  precession  is  small,  (ii)  the  viscosity  of  the 
liquid  is  small  (the  Reynolds  number  is  large),  ( i i i )  the  ratio  of  the  mass  of 
the  liquid  to  the  mass  of  the  solid  parts  of  the  system  is  small,  and  (iv)  the 


1 .  Engineering  Design  Handbook.  Liquid-Filled  Projectile  Design.  AMC  Pamphlet 
No.  706-165 ,  U.S .  Army  Materiel  Development  and  Readineee  Cornand, 
Washington ,  DC,  April  1969.  AD  853719 

2.  Greenspan,  H.P.,  The  Theory  of  Rotating  Fluids .  Cambridge  Univ.  Press,  NY, 
1969,  pp.  81  ff. 

3.  Stewartson,  K.,  "On  the  Stability  of  a  Spinning  Top  Containing  Liquid", 

J.  Fluid  Mesh..  Vol .  5,  No.  4,  September  1959,  pp.  577-592. 

4.  Phillips,  O.M.,  "Centrifugal  Waves",  J .  Fluid  Mech. .  Vol.  7,  1960,  pp. 
340-352. 

5.  Wedemeyer,  E.H.,  " Viscous  Corrections  to  Stewartson' s  Stability 
Criterion",  U.S.  Army  Ballistic  Research  Laboratory/ARRADCOM  Report 
No. 1325,  Aberdeen  Proving  Ground,  Maryland,  June  1966.  AD  489687 . 

6.  Wedemeyer,  E.H.,  "Viscous  Corrections  to  Stewartson' s  Stability 
Criterion" ,  AG ARP  Conference  Proceedings  No.  10.  NATO  Advisory  Group  for 
Aerospace  Research  and  Development,  Paris,  pp.  99-116,  September  1966. 


the  liquid  is  in  a  basic  state  of  riqid  rotation,  the  only  departure  being  the 
small  motion  due  to  the  inertial  oscillations.  In  addition,  the  influence  of 
the  liquid  on  the  system  frequency  and  qrowth  rate  is  approximated  as  being 
entirely  due  to  one  of  the  inertial  modes,  which  requires  that  the  system 
frequency  be  very  near  the  eigenfrequency  of  that  mode,  i.e. ,  that  the  system 
be  very  near  resonance. 

Gerber7  has  produced  a  linear  theory  in  which  the  contributions  of 
several  eigenmodes  are  summed,  and  which  is  therefore  valid  away  from 
resonance.  He  also  makes  no  assumptions  about  the  relative  masses  of  the 
liquid  and  solid  parts  of  the  system.  The  assumptions  of  small  precession 
angle,  small  viscosity,  and  a  basic  state  of  rigid  rotation  of  the  liquid 
remain,  and  there  is  the  important  new  requirement  that  the  liquid  fill  the 
cavity  completely.  For  those  cases  which  have  been  examined,  Gerber's  pre¬ 
dicted  growth  rates  agree  with  the  Stewartson-Wedemeyer  theory  in  the  regions 
where  both  theories  should  be  valid7. 

The  author  is  unaware  of  any  theory  predicting  the  motion  of  the  system 
when  the  liquid  is  not  in  a  basic  state  of  rigid  rotation,  for  example,  when 
it  is  being  spun  up.  However,  two  theories  have  been  presented  predicting  the 
eigenmodes  of  oscillation  for  this  case.  Lynn8  has  calculated  eigenfre- 
quencies  of  axisymmetric  (i.e.,  azimuthal ly-constant)  inertial  modes,  using  an 
inviscid  formulation.  Kitchens,  Gerber,  and  Sedney9  have  calculated  eigen¬ 
values  for  both  axisymmetric  and  nonaxi symmetric  modes,  with  a  viscous  formu¬ 
lation  valid  for  large  Reynolds  number. 

All  of  the  theories  mentioned  above  are  linear  in  the  angle  of  preces¬ 
sion,  e.  Scott10  used  a  formal  expansion  of  the  cavity  length  in  powers 
of  e  to  obtain  a  first-order  finite  amplitude  correction  to  the  inviscid 
eigenfrequencies.  Although  the  applicability  of  his  model  to  real  fluids 
seems  somewhat  doubtful,  since  he  included  no  viscous  effects,  his  predicted 


7.  Gerber,  N.,  "Moment  on  a  Liquid-Filled  Projectile:  Solid-Body  Rotation", 
U.S.  Army  Ballistic  Research  Laboratory /ARRADCOM  Report  (to  be  published) • 

8.  Lynn,  Y.M.,  "Free  Oscillations  of  a  Liquid  During  Spin-Up ",  U.S.  Army 
Ballistic  Research  Laboratory /ARRADCOM  Report  No.  1663 ,  Aberdeen  Proving 
Ground,  Maryland,  August  1973.  AD  769710. 

9.  Kitchens,  Jr.,  C.W.,  Gerber,  N.,  and  Sedney,  R.,  " Oscillations  of  a  Liquid 
in  a  Rotating  Cylinder:  Part  I.  Solid  Body  Rotation ",  U.S.  Amy 
Ballistic  Research  Laboratory /ARRADCOM  Report  ARBRL-TR-02081 ,  Aberdeen 
Proving  Ground,  Maryland,  June  1978.  AD  A0S77S9. 

10.  Scott,  W.E.,  "The  Large  Amplitude  Motion  of  a  Liquid-Filled  Gyroscope  and 
the  Non-interaction  of  Inertial  and  Roseby  Waves ",  J.  Fluid  Manh.s  Vol. 
72,  Part  4,  pp.  649-660,  1975. 
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frequency  corrections  did  seem  to  match  the  observed  shift  in  system  frequency 
at  maximum  growth  rate  observed  by  Scott  and  D'Amico11. 

Relevant  experimental  investigations  fall  into  two  broad  classes:  (1) 
studies  of  the  behavior  of  the  solid-liquid  systems;  and  (2)  studies  of  the 
motion  of  the  liquid  only.  Ward12  conducted  the  first  reported  experiment  of 
the  first  class,  in  which  he  crudely  measured  the  growth  rate  of  the  preces¬ 
sion  angle  of  a  liquid-filled  gyrostat,  finding  one  of  Stewartson's  predicted 
resonances.  Karpov1  conducted  a  series  of  experiments  using  liquid- fil led 
projectiles,  demonstrating  the  existence  of  instabilities  and  showing  quali¬ 
tatively  the  effects  of  viscosity  and  of  the  ratio  of  liquid  mass  to  solid 
mass.  Karpov  and  others  conducted  gyrostat  experiments,  recently  reviewed  by 
Whiting  and  Gerber13,  the  results  of  which,  for  a  very  limited  set  of  condi¬ 
tions,  agreed  moderately  well  with  the  predictions  of  the  Stewartson-Wedemeyer 
theory.  In  all  of  these  experiments  the  liquid  was  in  a  basic  state  of  rigid 
rotation. 

Karpov1  and  more  recently  Mark1*4  have  found  evidence  in  tests  with 
liquid-filled  projectiles  strongly  suggesting  that  system  instabilities  can 
occur  when  the  liquid  payload  is  being  spun  up.  D'Amico  and  Miller15  have 
observed  an  instability  of  a  projectile  filled  with  a  highly-viscous  liquid 
(Reynolds  number  about  1);  it  seems  unlikely  that  the  mechanism  of  resonance 
of  inertial  oscillations  with  projectile  precession  had  anything  to  do  with 
this  instability,  but  the  observations  make  the  practical  point  that  instabil¬ 
ities  at  smaller  viscosities  cannot  necessarily  be  made  to  go  away  by  simply 
increasing  the  viscosity. 

In  all  of  the  experiments  using  liquid-filled  gyrostats,  it  was  found 
that  the  growth  rates  of  instabilities  deviated  from  linear  dependence  on  t  at 
rather  small  angles.  Scott  and  D'Amico11  found,  for  example,  that  growth 


11.  Scott ,  W.E.,  and  D'Amico,  W.P.,  "Amplitude-Dependent  Behavior  of  a 
Liquid-Filled  Gyroscope",  J.  Fluid  Mech.,  Vol.  60,  Part  4,  pp.  751-758 , 
1973. 

12.  Ward,  B.N.,  Appendix  to  Reference  3. 

13.  Whiting,  R.D.,  and  Gerber,  N.,  "Dynamics  of  a  Liquid-Filled  Gyroscope: 
Update  of  Theory  and  Experiment",  U.S.  Army  Ballistic  Research 
Laboratory /ARRADCOM  Technical  Report  ARBRL-TR-02221 ,  Aberdeen  Proving 
Ground,  Maryland,  March  1980.  AD  A083886. 

14.  Mark,  A.,  "Measurements  of  Angular  Momentum  Transfer  in  Liquid-Filled 
Projectiles",  U.S.  Army  Ballistic  Research  Laboratory /ARRADCOM  Report 
ARBRL-TR-2029,  Aberdeen  Proving  Ground,  Maryland,  November  1977.  AD 
A051342. 

15.  D'Amico,  W.P.,  and  Miller,  M.,  "Flight  Instability  Produced  by  a  Rapidly- 
Spinning,  Highly-V iscoue  Liquid",  «/.  Spacecraft,  and  Rncka+.*f  Vol.  16,  No. 
1,  pp.  62-64,  January  1979. 
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rates  of  instabilities,  and  the  system  frequency  at  which  maximum  growth  rate 
occurred,  changed  at  precession  anqles  as  small  as  9  *  10" 3  radians. 

The  first  quantitative  investigation  of  the  second  class  (studies  of  the 
motion  of  the  liquid  only)  was  by  Fult218,  who  excited  axisymmetric  modes  in 
a  cylinder  by  means  of  axial  oscillations  of  a  disk  in  the  liquid.  He 
measured  the  frequencies  of  forcing  which  produced  the  greatest  response,  as 
observed  by  means  of  dye  injected  into  the  liquid.  For  Reynolds  numbers 
around  5  *  103,  he  found  agreement  to  within  one  percent  with  eigenvalues 
computed  using  inviscid  theory. 

Thompson17  used  gravitational  forcing  of  the  free  surface  of  the  liquid 
in  a  tilted  cylinder  to  drive  a  nonsymmetric  mode  of  frequency  zero  (in  the 
laboratory  frame  of  reference).  He  found  that  nonlinear  effects  produced  an 
azimuthal ly-constant,  radially-varying  secondary  flow,  which  could  suffer  a 
shear  instability  if  the  liquid  depth  were  such  that  an  inertial  mode  had 
natural  frequency  near  zero. 

McEwan18  forced  nonsymmetric  oscillations  by  causing  the  rotation  axis  of 
a  tilted,  rigid  lid  to  cone  about  the  rotation  axis  of  the  cylinder,  at  a 
speed  different  from  that  of  the  cylinder  rotation.  Using  flow  visualization, 
he  demonstrated  the  presence  of  a  number  of  nonsymmetric  modes  near  the  criti¬ 
cal  aspect  ratios  predicted  by  linear,  inviscid  theory.  He  was  not  able  to 
show  viscous  effects  on  the  natural  modes  directly.  He  did  show  the  presence 
of  a  mean  secondary  flow  apparently  qualitatively  similar  to  that  seen  by 
Thompson.  Although  he  did  not  resolve  the  radial  structure  of  this  flow,  some 
of  his  measurements  roughly  indicated  the  change,  caused  by  the  secondary 
flow,  in  the  critical  cavity  dimensions  required  for  resonance  of  a  natural 
mode. 


Both  Thompson  and  McEwan  observed  flow  instabilities  whenever  the  distur¬ 
bance  angle  (in  McEwan' s  case  the  tilt  of  the  lid,  in  Thompson's  case  the  tilt 
of  the  cylinder)  was  sufficiently  large.  Thompson  found  that  the  flow  became 
unstable  at  a  disturbance  angle  of  about  10”3  radians,  when  the  fundamental 
mode  was  in  resonance,  and  about  0.015  radians  when  the  next  higher  mode  was 
in  resonance.  The  only  value  of  Reynolds  number  he  used  was  6.6  *  104. 

McEwan  found  flow  instability  at  a  tilt  angle  of  .018  radians,  for 
Reynolds  number  greater  than  9  «  104.  He  presents  data  indicating  that  higher 
viscosity  inhibits  the  instability.  He  remarked  that  "Collapse  [instability] 


16.  Fultz,  D.,  "A  Note  on  Overstability  and  the  Elastoid- Inertia  Oscillations 

of  Kelvin,  Solberg,  and  Bjerkness",  J.  Meteorology,  Vol.  16,  pp.  199-208, 
April  1969.  ~~  ~~ 

17.  Thompson,  R.,  "Diurnal  Tides  and  Shear  Instabilities  in  a  Rotating 
Cylinder",  J.  Fluid  Mech.,  Vol.  40,  Part  4,  pp.  737-751,  1970. 

18.  McEwan,  A.D.,  '"Inertial  Oscillations  in  a  Rotating  Fluid  Cylinder”,  J. 

Fluid  Mech.,  Vol.  40,  Part  3,  pp.  603-640,  1970.  ~~ 
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would  also  be  avoided  in  all  other  configurations  for  sufficiently  weak 
excitation  rate  or  low  rotation  rate  [low  Reynolds  number]...",  but  he  did  not 
report  what  angles  or  rotation  rates  were  "sufficiently  weak". 


Aldridge19’20  forced  axisymmetric  oscillations  in  a  cylinder  during  spin- 
up  of  the  liquid  by  causing  the  rotational  speed  of  the  cylinder  to  vary 
periodically.  He  determined  the  response  of  the  liquid  by  measuring  the 
fluctuating  pressure  in  the  cylinder,  and  found  that  as  time  progressed,  that 
is,  as  the  liquid  spun-up,  the  amplitude  of  the  pressure  fluctuations  came  to 
a  peak  value,  then  decayed.  The  time  at  which  the  peak  occurred  was  usually 
later  than  the  time  at  which  Lynn8  predicted  that  the  liquid  eigenfrequency 
would  match  the  forcing  frequency,  being  about  fifteen  percent  later  at  most. 
Subsequently,  Gerber  and  Sedney  (personal  communication)  have  used  the 
solution  of  Reference  9  to  make  predictions  analogous  to  those  of  Lynn.  Their 
results  lie  between  Lynn's  predictions  and  Aldridge's  data.  Stergiopoulos  and 
Aldridge21  have  studied  nonsymmetric  oscillations  during  spin-up  of  a  liquid 
in  a  slightly  tilted  cylinder  with  the  upper  end  open,  the  same  type  set-up 
used  by  Thompson  17.  They  were  able  to  determine  the  time  evolution  of  reso¬ 
nant  height  as  the  liquid  progressed  to  rigid  rotation,  a  result  not  yet 
available  in  theory.  When  the  liquid  was  in  rigid  rotation,  they  found  non¬ 
linear  effects  at  the  smallest  tilt  angle  studied  (.005  radians),  and  insta¬ 
bility  of  the  flow  at  any  angle  above  .01  radians,  when  the  depth  of  the 
liquid  was  such  that  the  second  lowest  mode  was  in  resonance  (the  same  mode 
for  which  Thompson  found  instability  above  .015  radians). 

It  is  useful  to  summarize  the  "state-of-the-art"  as  represented  by  the 
theories  and  experiments  described  above.  In  regard  to  the  motion  of  the 
liquid,  it  can  be  said  that  the  existence  of  discrete  frequencies  of  oscil¬ 
lation  in  the  liquid1’2’9  is  amply  confirmed16"18.  For  liquids  of  small 
viscosity,  the  spectrum  of  allowed  frequencies  seems  to  have  nearly  the 
dependence  on  geometry  predicted  by  linear,  inviscid  theory,  for  the  lower 
wave  number  modes  of  oscillation16  19’21.  The  effects  of  viscosity5’6’9  on 
the  oscillations  have  not  been  quantitatively  confirmed,  though  it  is  clear 
that  there  are  effects  18 ,  and  that  in  some  respects  they  scale  with  viscosity 
more-or-less  as  viscous  theory  predicts.  The  amplitude  attained  by  forced 
nonsymmetric  oscillations  has  not  been  accurately  measured,  except  for  the 
mode  with  frequency  zero21.  The  existence  of  nonlinear  effects  on  the 


19.  Aldridge,  K.,  "Experimental  Verification  of  the  Inertial  Oscillations  of 
a  Fluid  in  a  Cylinder  During  Spinr-Up",  U.S.  Army  Ballistic  Research 
Laboratory /ARRADCOM  Contract  Report  No.  273,  Aberdeen  Proving  Ground, 
Maryland,  November  1975.  AD  #A018797 

20.  Aldridge,  K.,  "Axisymmetric  Inertial  Oscillations  of  a  Fluid  in  a 
Cylindrical  Cavity  During  Spin-Up  from  Rest",  Geophysical  and 
Ast.mphyfsi.ruil  Fluid  Dyrumi.r.s,  Vol.  8,  No.  4,  pp.  279-301,  1977. 

21.  Stergiopoulos,  S.,  and  Aldridge,  K.,  "Inertial  Waves  in  a  Partially- 
Filled  Cylindrical  Container  during  Spin-Up  from  Rest",  submitted  to 

/T/i/MnluMx'  /tm/7  A  nmt/rnvt/tn 


response  of  a  nonsymmetric  mode  and  a  change  in  the  mean  motion  of  the  liquid 
due  to  these  nonlinearities  have  been  well  demonstrated  at  angles  as  small  as 
.001  radians17’18.  It  has  been  established  that  the  flow  may  become  unstable, 
possibly  due  to  shear  instability  of  the  secondary  flow,  at  angles  of  the 
order  of  .001  radians17’ 18._  It  has  been  shown  that  inertial  oscillations  can 
be  driven  during  spin-up19’21,  and  that  for  several  axi symmetric  modes,  and 
some  nonsymmetric  modes  of  zero  frequency,  the  trajectories  of  eigen frequency 
in  time  are  smooth.  For  the  nonsymmetric  modes,  nonlinear  effects  have  been 
found  during  spin-up21.  For  the  axisymmetric  modes,  the  predicted 
trajectories  of  eigenfrequency  are  reasonably  close  to  the  measured 
trajectories  of  peak  response19’20,  though  no  demonstration  has  been  made  that 
the  eigenfrequency  calculated  by  the  theories8’9  should  be  identified  with  the 
peak  response. 

With  respect  to  system  response,  there  has  been  only  one  complete  check 
of  the  predicted  versus  measured  precession  frequency  and  growth  rate  of  a 
liquid- filled  gyrostat13.  Several  partial  checks  are  available,  including 
cases  where  the  liquid  does  not  completely  fill  the  cavity.  These  experiments 
match  the  theoretical  predictions  to  within  ten  to  twenty-five  percent13. 

There  is  clear  evidence  in  gyrostats  of  an  effect  of  finite  amplitude  of 
precession  angle  on  the  growth  rate  at  angles  as  small  as  .009  radians11,  but 
it  cannot  be  said  whether  the  effect  is  due  primarily  to  nonlinear  processes 
in  the  liquid  or  in  the  gyrostat  response.  The  point  that  liquid  motions  seem 
to  become  unstable17’ 18  at  smaller  angles  than  those  at  which  the  growth  rates 
of  gyrostats  change  has  not  been  addressed.  System  instabilities  can  occur 
while  the  liquid  is  being  spun-up1'14,  but  such  instabilities  cannot  be 
predicted. 


B.  Purpose 

The  objective  of  this  experiment  was  to  measure  the  flow  due  to  forced 
nonsymmetric  oscillations  when  the  liquid  was  in  a  basic  state  of  rigid 
rotation,  by  measuring  pressure  fluctuations,  using  a  modified  version  of 
Aldridge's  technique.  The  effects  of  viscosity,  fill  ratio,  and  geometry  were 
examined,  for  a  range  of  coning  frequencies  in  the  region  applicable  to 
liquid-filled  projectiles.  For  each  combination  of  these  parameters,  the 
limiting  precession  angle  at  which  linear  theory  could  be  used  was  determined, 
as  were  the  nature  of  the  nonlinear  effects  when  this  angle  was  exceeded. 

In  addition,  apparatus  and  techniques  were  developed  which  could  be  used 
to  conduct  similar  studies  of  the  flow  while  the  liquid  is  being  spun-up  from 
rest,  these  studies  to  form  a  separate,  follow-on  effort. 


II.  EXPERIMENTAL  METHOD 

A.  Diagnostic  Variable 

Fluctuating  pressure  was  chosen  as  the  quantity  to  be  measured.  It  is  a 
primitive  flow  variable,  and  is  the  quantity  whose  value  is  required  for 
computing  overturning  moments  on  a  container.  It  was  decided  to  measure  the 
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pressure  at  an  end  wall,  so  that  the  measuring  device  need  not  intrude  into 
the  flow.  The  inference  of  interior  pressure  from  measurements  at  a  wall 
requires  that  the  boundary  layer  approximation  be  valid,  so  that  the  normal 
derivative  of  pressure  near  the  wall  is  small.  This  in  turn  requires  that  the 
Reynolds  number  be  large. 

The  choice  of  an  end  wall  in  preference  to  the  side  wall  was  based  purely 
on  convenience.  One  might  argue  that  the  pressure  on  the  side  wall  contri¬ 
butes  more  to  the  overturning  moment  exerted  on  the  container  and  is  therefore 
more  relevant  than  pressure  on  the  end  wall.  However,  the  pressure  measured 
anywhere  in  the  liquid  should  be  equally  useful  as  a  flow  diagnostic.  The 
only  new  information  to  be  gained  from  sidewall  pressure  is  spatial  structure, 
which  is,  of  course,  important.  If  it  is  later  decided  that  measurements  at 
the  sidewall  are  desirable,  they  can  be  made  with  minor  changes  in  the  appar¬ 
atus. 

B.  Independent  Variables 

The  quantities  varied  during  the  experiment,  and  the  apparatus  to  imple¬ 
ment  the  variations,  were: 

1.  Cylinder  Geometry.  The  cylinder  (Figure  1)  was  constructed  of 
Lucite,  with  a  mean  inner  radius,  a,  of  31.761+  .012mm.  Aluminum  top  and 
bottom  inserts  of  several  different  lengths  were  made,  so  that  the  interior 
length,  2c,  could  be  changed  by  changing  inserts.  Two  lengths  were  used  in 
this  experiment:  199.972+  .010mm,  and  66.756+  .011mm.  The  aspect  ratios, 
c/a,  were  thus  3.1481  and  1.0509  at  22°  C.  The  total  uncertainty  in  c/a, 
including  the  effect  of  temperature  fluctuations,  is  estimated  to  be  .0012  for 
c/a  =  3.1481,  and  .0004  for  c/a  *  1.0509. 

2.  Coning  Angle.  The  cylinder  was  mounted  in  a  gimbal  case  (Figure  2) 
which  contained  a  motor  to  maintain  the  rotational  speed  ft.  The  case  had  a 
shaft  extending  from  its  lower  end,  which  mated  with  a  bushing  gripped  by  an 
adjustable  cam  (Figure  3).  The  cam  was  rotated  by  a  second  electric  motor;  as 
it  rotated,  it  caused  the  cage,  hence  the  axis  of  rotation  of  the  cylinder,  to 
cone  about  a  line  passing  through  the  center  of  the  cam  and  the  pivot  point  of 
the  cylinder  (the  "coning  axis"  in  Figures  2  and  3).  The  pivot  point,  that 
is,  the  intersection  of  the  axes  of  the  inner  and  outer  gimbal  bearings,  was 
coincident  with  the  center  of  the  cylinder.  The  coning  angle  e  was  adjusted 
by  adjusting  the  offset  screw  in  the  cam.  The  angle  was  measured  (after  the 
locking  screw  was  tightened)  in  one  of  two  ways:  (1)  if  the  coning  angle  was 
less  than  3.5  x  10"3  rad.,  the  displacement  of  a  point  on  the  motor  mount  was 
measured_  with  a  dial  indicator;  (2)  if  the  coning  angle  was  greater  than 
3.5  x  10" 3  rad.,  the  number  of  turns  of  the  offset  screw  used  to  set  the  angle 
was  recorded.  Either  measurement  could  be  turned  into  horizontal  displacement 
of  a  known  point  on  the  gimbal  cage;  with  the  previously-measured  distance  of 
that  point  from  the  pivot  point,  e  could  then  be  determined.  Geometric 
relations_ accurate  to  first  order  in  e  were  used,  so  that  e  was  measured  to 
+1.3  x  10  4  rad.  when  it  was  less  than  3.5  x  10"3  rad.,  or  +9.7  x  10"6  rad., 
when  it  was  greater  than  3.5  x  10" 3  rad. 

The  smallest  coning  angle  used  was  about  1.7  x  10"4  rad.;  the  largest  was 
about  3.5  x  io"2  rad. 
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3.  Viscosity  Three  different  liquids  were  used:  Dow  Corning  200 

Silicon  Fluids  with  nominal  kinematic  viscosities  v  of  1,  5,  and  100 

centistokes.  Details  of  the  variation  of  viscosity  and  density  of  the  liquids 
with  temperature  are  described  in  Appendix  A.  Two  cylinder  rotational  speeds 
were  used:  523.6  rad./s,  and  433.3  rad./s*.  Two  non-dimensional  numbers 
involving  v  are  commonly  used:  the  Reynolds  number.  Re  =  i 2a2/v,  and  the  Ekman 
number,  E  5  v/flc2.  In  this  experiment,  Re  was  varied  from  about  5  x  io3  to 
5  x  10s;  E  varied  from  about  2  x  io"7  to  2  *  10" b. 

4.  Fi  11  Ratio  In  addition  to  the  filled  case,  data  were  taken  with 
about  50  m'l  {about  eight  percent  of  the  cavity  volume)  of  liquid  removed  from 
the  cylinder  with  aspect  ratio  3.1481.  This  gave  a  fill  ratio,  f  {=  liquid 
volume/cavity  volume)  of  about  0.92.  If  the  axial  acceleration  of  the  liquid 
(e.g.,  gravity)  were  zero,  then  the  interface  between  the  liquid  and  the  air 
would  be  cylindrical  (neglecting  distortions  due  to  oscillations).  In  this 

case,  the  air  core  radius,  b,  would  be  a(l-f)^,  or  about  8.9mm.  The  appro¬ 
priate  measure  of  the  axial  acceleration  is  the  Froude  number,  F  (h  U2b/g) , 
which  was  about  172  for  this  experiment.  This  caused  the  air-liquid  interface 
to  be  a  truncated  paraboloid,  with  the  radius  at  the  upper  and  lower  ends 
respectively  6.3%  greater  and  6.7%  less  than  the  value  of  b  computed  above.  F 
was  not  varied  during  the  experiment. 

5.  History  of  Cylinder  Spin.  In  this  experiment,  u  was  held  constant  to 
within  about  one  percent  for  at  least  five  minutes  before  data  were  recorded, 
ensuring  that  the  liquid  had  come  to  a  basic  state  of  rigid  rotation.  While 
data  were  actually  being  recorded,  a  was  kept  constant  to  within  .05%. 

6.  History  of  Coning  Frequency.  The  coning  frequency  was  kept  constant 
to  within  about  0.1%  long  enough  for  transients  to  die  away  (see  paragraph  D6 
below). 

C.  Instrumentation 


The  principal  instrumentation  included  the  following: 

1.  Pressure  Transducers.  Kulite  XT-39-190-25A  absolute  transducers  were 
used.  These  "transducers  were  specially  built  with  metal  stops  behind  the 
sensing  diaphram.  Large  deflections  which  would  ordinarily  damage  the 
transducers  are  not  permitted  due  to  the  metal  stop.  Over  pressure  ratios  of 
40/1  are  permitted  without  a  loss  in  calibration.  Two  transducers  were  flush 
mounted  in  the  bottom  end  face  of  the  cavity,  as  shown  in  Figure  4.  One 
transducer  was  at  a  non-dimensional  radius  r/a  =  0.436,  the  other  at  r/a  = 
0.668;  the  transducers  differed  in  azimuthal  location  by  about  u/2.  The 
sensing  area  of  each  transducer  had  a  diameter  of  about  2.16mm,  or  0.068a. 
The  transducer  signals  were  fed  into  amplifiers 


*The  lower  speed  was  used  in  the  cases  of  partial  fill,  to  avoid  a  resonance 
with  the  supporting  structure. 
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which  provided  a  gain  of  about  700  in  the  range  70-90  Hz  (the  range  of  concern 
in  this  experiment);  the  amplifier  cut  off  frequencies  were  slightly  below  10 
Hz.  Details  of  the  calibration  of  the  transducers  and  amplifiers  are  in 
Appendix  B. 

2.  Telemetry.  The  pressure  signals  were  transferred  from  the  rotating 
cylinder  frame  of  reference  to  the  laboratory  frame  by  a  telemetry  link,  shown 
schematically  in  Figure  5.  The  outputs  of  the  transducer  amplifiers  were 
input  to  two  Omnitek  Model  22A1E  voltage  controlled  oscillators  (V COs)  with  40 
kHz  and  52.5  kHz  center  frequencies,  whose  outputs  were  qomblned  and  sent  to 
an  Optimax  Model  4001  FM  transmitter.  The  transmitted  signal  was  received, 
amplified  by  a  Mu-Del  Electronics,  Inc.,  Model  MDA-2132  RF  preamplifier,  and 
sent  to  a  Defense  Electronics,  Inc.,  Model  TMR-511  receiver.  The  video  output 
from  the  receiver  passed  to  two  EMR  Model  4142  tunable  discriminators,  set  for 
40  kHz  and  52.5kHz,  whose  outputs  were  then  equal  to  the  VCO  inputs  (less  an 
adjustable  offset).  The  overall  transfer  function  of  all  of  this  equipment, 
including  the  transducer  amplifiers,  was  determined  using  a  Hewlett-Packard  HP 
3582A  spectrum  analyzer.  Details  of  the  equipment  settings  and  calibration 
procedures  are  shown  in  Appendix  B. 

3.  Signal  Analyzer.  The  same  HP  3582A  spectrum  analyzer  was  used  to 
analyze  the  pressure  signals.  This  device  computed  discrete  Fourier  trans¬ 
forms  of  two  512-point  time  samples  of  data  gathered  simultaneously,  enabling 
data  to  be  taken  on  two  channels  at  once.  The  instrument  was  set  to  take 
evenly-spaced  samples  over  a  12.5-second  interval;  low-pass  filters  within  the 
instrument  prevented  aliasing.  The  instrument  performed  complex  multiplica¬ 
tions  of  the  signals  by  a  sinusoid  of  80  Hz  frequency,  which  enabled  the  128- 
point  discrete  transforms  to  cover  a  span  from  75  Hz  to  85  Hz.  Details  of  the 
procedure  are  in  the  manufacturer's  literature22. 

The  spectrum  analyzer  presented  both  magnitude  and  phase  of  the  trans¬ 
form,  and  both  were  recorded.  The  manufacturer  cites  accuracy  of  about  ±6%  in 
magnitude  and  +10°  In  phase  over  the  entire  operating  range  of  the  instrument; 
tests  with  controlled  sources  in  only  the  range  used  in  the  experiment  showed 
the  actual  worst  case  error  in  magnitude  to  be  about  1.5%.  The  Instrument 
offers  three  choices  of  data  windowing.  The  "flat  top"  window  was  selected, 
since  it  offers  the  best  amplitude  accuracy  for  signals  which  differ  in  fre¬ 
quency  from  any  Fourier  frequency22. 

The  spectra  obtained  during  the  experiments  had,  in  general,  two  dominant 
peaks,  as  shown  in  Figure  6.  One  was  at  the  spin  frequency,  and  was  only 
slightly  affected  in  its  amplitude  by  the  coning  frequency.  Its  amplitude  was 
reduced  by  a  factor  of  two  or  three  when  there  was  no  liquid  in  the  cylinder, 
and  was  increased  if  the  cylinder  was  poorly  balanced.  The  phases  of  the 
transforms  of  the  two  pressure  signals  at  this  frequency  were  equal,  indicat¬ 
ing  that  there  was  no  azimuthal  variation  of  this  component  in  the  rotating 
frame  of  reference.  This  peak  was  taken  to  be  due  to  disturbances  caused  by 


22.  " Operating  Manual,  Model  3S82A  Spectrum  Analyser",  HP  Manual  03582-90000, 

Hewlett  Packard  Co.,  Loveland,  Colorado,  1978. 
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vibration  due  to  residual  imbalance  and  bearing  imperfections,  plus  some 
electrical  interference  or  other  effect  due  to  azimuthal  variations  in  the 
pattern  of  the  transmitting  antenna. 

The  second  peak  was  at  a  frequency  equal  to  the  spin  frequency  minus  the 
coning  frequency,  which  is  what  is  expected  from  a  transducer  sweeping  at  the 
spin  frequency  through  a  wave  travelling  at  the  coning  frequency  (spin  and 
coning  were  always  in  the  same  directions).  This  peak  was  taken  to  be  the 
response  of  the  liquid  to  the  coning.  Its  amplitude  was  the  amplitude  datum 
recorded,  and  its  variation  with  coning  frequency  and  coning  angle  is  the 
subject  of  the  discussion  in  Section  III  of  this  report.  The  phases  of  the 
two  signals  at  this  frequency  differed  by  approximately  ir/2,  indicating  that 
there  was  azimuthal  variation  in  the  rotating  frame  of  reference.  Figure  7 
shows  the  amplitude  spectra  taken  on  one  channel  at  ten  different  coning 
frequencies,  overlaid  on  the  same  sole.  This  figure  clearly  illustrates  the 
large  variation  of  the  forced  response,  and  the  much  smaller  variation  of  the 
peak  at  the  spin  frequency. 

4.  Interval  Timers.  Magnetic  sensors  and  associated  signal  conditioners 
were  used  to  produce  electrical  pulses  on  three  circuits:  one  when  the 
cylinder  completed  one  revolution;  one  when  the  coning  cam  completed  one 
revolution  (hence  when  the  cylinder  axis  completed  one  cycle  of  coning);  and 
one  whenever  the  other  two  pulses  happened  to  occur  within  about  800  ps  of 
each  other.  This  last  pulse  thus  occurred  when  the  transducers  were  within 
.0175  radians  of  a  known  position  relative  to  the  plane  defined  by  the 
cylinder  rotation  axis  and  the  axis  about  which  it  was  coning,  that  is,  when 
the  phases  of  the  transducers  with  respect  to  the  coning  were  known.  This 
pulse  was  used  to  trigger  the  spectrum  analyzer  to  begin  gathering  a  sample  of 
data,  and  in  principle  allowed  the  computed  phase  of  the  pressure  signals  to 
be  related  to  the  phase  of  the  coning. 

The  pulses  produced  by  cylinder  rotation  were  fed  to  a  Tektronix  model 
DC503  counter  set  to  time  the  interval  between  pulses;  the  counter  displayed 
the  period  of  cylinder  rotation  to  the  nearest  microsecond.  This  signal  was 
observed  continuously  during  the  time  a  data  sample  was  being  gathered,  and 
the  constancy  of  cylinder  speed  cite,  earlier  comes  from  observations  of  the 
variation  of  the  rotation  period.  In  every  case,  the  actual  period  was 
recorded  to  the  nearest  10  us. 

The  pulses  produced  by  rotation  of  the  coning  cam  were  fed  to  the  BRL  Yaw 
Processor,  which  stored  each  period  of  revolution  in  a  digital  memory.  After 
each  sample  of  data  was  gathered  by  the  spectrum  analyzer,  a  HP  9830A  calcu¬ 
lator  recovered  the  stored  periods,  and  calculated  the  mean,  minimum,  and 
maximum  periods  for  the  interval  (typically  48  periods),  as  well  as  the  stand¬ 
ard  deviation  of  the  periods.  The  mean  period  was  recorded  to  the  nearest 
10  ps.  The  other  information  was  used  to  determine  the  extent  of  fluctuations 
in  coning  period,  which  was  used  to  compute  the  constancy  of  coning  frequency 
cited  previously. 

5.  Computer.  A  Hewlett-Packard  9845B  computer  was  connected  to  the 
spectrum  analyzer,  and  the  amplitudes  and  phases  of  the  two  pressure  signals 
at  the  frequency  of  the  forced  oscillation  were  transmitted  directly  from  the 
analyzer  to  the  computer.  Cylinder  rotation  period  and  mean  coning  period 
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were  entered  into  the  computer  by  hand.  The  computer  collected  all  of  the 
data  from  one  "run",  that  is,  one  value  of  all  parameters  except  coning 
period,  and  a  sequence  of  values  of  coning  period,  into  a  file.  All  subse¬ 
quent  manipulation  of  the  data  was  performed  on  this  computer. 

Although  it  was  not  done  in  this  experiment,  the  HP  9830A  calculator 
which  computed  the  coning  period  could  also  have  been  connected  to  the  HP 
9845B  computer.  The  addition  of  a  counter  with  the  appropriate  interface 
( I  EE-488)  and  two  digital-to-analog  converters  would  in  the  future  allow 
completely  automated  control  of  cylinder  and  coning  speed  plus  data  gathering. 

D.  Procedure 


The  procedure  for  one  run  was  as  follows: 

1.  The  cylinder  and  inserts  were  assembled,  and  the  cavity  was  filled 
with  liquid.  The  cavity  was  connected  to  a  vacuum  pump  for  (typically)  an 
hour,  to  deaerate  the  liquid.  Since  this  caused  substantial  cooling  of  the 
liquid,  the  apparatus  was  allowed  to  settle,  usually  overnight,  to  allow  the 
liquid  to  come  to  room  temperature. 

2.  The  coning  angle  was  set  and  measured,  as  discussed  previously. 

3.  The  electronic  instruments  were  turned  on  and  allowed  to  warm  up  for 
at  least  30  minutes.  The  discriminators  were  adjusted  to  obtain  unity  gain  in 
the  telemetry  link  (the  procedure  is  detailed  in  Appendix  B). 

4.  Room  temperature  was  measured  and  recorded.  It  was  assumed  that  the 
liquid  thereafter  remained  at  room  temperature,  and  this  value  was  used  later 
to  determine  the  density  and  viscosity  of  the  liquid. 

5.  The  cylinder  was  inserted  in  the  gimbal  cage,  and  brought  up  to 
speed.  At  least  five  minutes  elapsed  before  any  data  were  taken.  Since  the 

1  /2 

longest  linear  spin-up  time  (c/(siv)  '  )  for  this  experiment  was  4.8  seconds, 
this  delay  ensured  that  the  liquid  was  fully  spun  up. 

6.  The  coning  cam  was  brought  to  its  selected  speed.  The  spectral 

display  was  monitored  to  determine  when  the  amplitude  of  the  forced  oscil¬ 
lations  had  stabilized,  or  for  at  least  30  seconds.  In  every  case  in  which 
the  display  did  stabilize,  that  occurred  within  30  seconds:  therefore  in 
practice  the  time  limit  governed. 

7.  If  the  display  stabilized,  a  data  sample  was  taken.  In  many  earlier 
cases  in  which  the  display  did  not  stabilize  (the  significance  of  which  is 
discussed  in  Section  III),  several  data  samples  were  taken.  In  later  cases, 
if  the  display  did  not  stabilize,  the  spectrum  analyzer  was  switched  to  a  2.5- 
second  sampling  time,  and  the  maximum  and  minimum  values  attained  over  some 
period  of  time,  typically  30  seconds,  were  recorded. 

Steps  6  and  7  were  repeated  for  all  selected  coning  frequencies. 
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Because  the  data  Mere  immediately  available  in  the  computer,  plots  of 
suitably  non-dimensional i zed  quantities  could  be  vieMed  within  a  few  minutes 
of  the  completion  of  a  run,  either  singly  or  overlaid  with  plots  from  previous 
runs.  This  facility  proved  invaluable  in  guiding  the  choice  of  conditions  of 
future  runs. 

III.  RESULTS 

A.  Comments  on  the  Measurements 

1.  Measurement  of  Phase.  The  effort  to  measure  phase  was  not  success¬ 
ful.  For  reasons  associated  with  the  characteristics  of  the  data  window 
employed  by  the  spectrum  analyzer,  errors  in  the  measurement  of  ft  and  w  caused 
much  larger  relative  errors  in  the  measured  phase.  The  estimated  errors  of 
.05%  in  ft  and  .1%  in  u>  can  be  shown  to  produce  errors  of  about  «/2  radians  in 
phase,  and  the  scatter  of  the  actual  data,  as  much  as  w  radians,  more  typi¬ 
cally  2w/3  radians,  in  nominally  identical  cases,  is  consistent  with  the  cited 
errors  in  ft  and  w. 

2.  Measurements  of  Amplitude.  The  results  of  the  pressure  measurements 
were  much  better  than  had  been  expected.  This  proved  extremely  fortunate,  as 
it  became  necessary  to  make  measurements  at  much  smaller  coning  angles  than 
had  been  anticipated  when  the  equipment  was  designed. 

The  results  of  the  measurements  will  be  presented  in  the  form  of  a  dimen¬ 
sionless  pressure  coefficient, 

A 

cp  =  p/epu2a2, 

A 

where  p  is  the  amplitude  of  the  fluctuating  part  of  the  pressure:  p  = 

A 

constant  +  p  sin  (wftt  +  <t>).  The  accuracy  with  which  Cp  could  be  determined 

depended  on  the  level  of  the  signal  and  on  the  coning  angle  e.  Using  the 
values  of  Cp  presented  later,  the  previously-cited  uncertainties  in  e,  ft,  a, 

and  the  electronic  signals,  the  relative  uncertainty  in  Cp  was  determined,  and 

is  presented  in  Figure  8  for  several  values  of  e.  For  a  few  points  in  the 
plots  of  Cp  presented  in  the  following  sections,  error  bars  are  indicated. 

These  reflect  the  relationship  shown  in  Figure  8. 

B.  General  Results 

For  the  purpose  of  discussing  the  data,  four  definitions  are  helpful: 

1.  "Linear  behavior":  linear  variation  of  p  with  e,  so  that  Cp  does  not 
depend  on  e; 

A 

2.  "Nonlinear  behavior":  nonlinear  variation  of  p  with  e,  so  that  Cp 
varies  with  e; 

A 

3.  "Nonstationary  behavior":  variation  of  p  with  time,  hence  of  Cp  with 
time  (thus  included  in  nonlinear  behavior);  and 


linear  regime".  In  which 


4.  "Regime":  a  part  of  parameter  space,  e.g.,  " 
some  specified  behavior  occurs. 

The  major  result  of  this  experiment  is  that  linear,  nonlinear,  and  non¬ 
stationary  regimes  were  all  seen.  This  is  illustrated  in  Figure  9,  which 
presents  Cp  plotted  against  <*>  for  one  choice  of  aspect  ratio,  Reynolds  number, 

and  fill  ratio.  The  different  symbols  represent  data  taken  at  different 
coning  angles.  Symbols  connected  by  a  vertical  line  represent  measurements 
for  which  the  amplitude  of  the  spectral  peak  of  the  forced  oscillation  did  not 
settle  to  a  constant  value,  i.e.,  nonstationary  behavior.  The  locations  of 
the  symbols  at  the  ends  of  the  line  are  Cp  computed  from  the  extreme  values  of 

the  spectral  peak.  They  must  not  be  interpreted  as  extremes  of  Cp,  however, 

since  the  spectrum  analyzer  operates  on  a  finite-length  sample,  thus  averaging 
over  the  sample  time  (up  to  12.5  seconds,  depending  on  the  frequency  span)  and 

since  no  particular  effort  was  made  in  every  case  to  follow  the  fluctuations 

in  time.  The  proper  interpretation  is  that  these  points  simply  indicate 
nonstationary  behavior. 

Where  different  symbols  lie  along  i**  same  curve,  linear  be_havior  is 
indicated.  In  Figure  9,  this  can  be  seen  to  occur  for  e  <  3.51  x  io“4.  For 
e  =  7.78  *  10~4,  Cp  has  a  constant  value  (*n  time),  but  near  its  peak  it 

diverges  from  Cp  at  smaller  e.  This  indicates  nonlinear  but  not  nonstationary 

behavior. 

In  every  case  studied  in  this  experiment,  Cp  attained  a  peak  value  for 

some  u.  Comparisons  of  the  w' s  of  peak  Cp  with  their  theoretical  values  for 

the  various  cases  are  presented  later,  but  as  a  general  comment,  there  was 
never  sufficient  disagreement  between  measured  and  theoretical  values  to  cloud 
the  identity,  that  is,  the  mode  numbers,  of  the  dominant  inertial  mode  near 
the  peak. 

C.  The  Linear  Regime 

All  of  the  theories  mentioned  in  Section  I  of  this  report,  except  that  of 
Scott10,  are  linear  in  e.  The  question  of  the  of  the  theories  is 

properly  addressed  by  comparing  their  predictions  with  data  which  show  linear 
behavior.  Nonlinear  behavior  then  provides  the  limits  of  .In 

this  section,  except  where  otherwise  noted,  only  data  in  the  linear  regime  are 
presented.  Differences  between  points  at  the  same  aspect  ratio,  Reynolds 
number,  and  fill  ratio  are  felt  to  represent  scatter. 

Representative  error  estimates  are  indicated  for  a  few  points.  In  most 
cases,  the  apparent  scatter  of  the  data  is  consistent  with  the  error 
estimates.  For  the  case  c/a  =  1.0509,  Re  =  5  x  105,  f  =  1.00,  shown  in  Figure 
10(d),  the  scatter  between  runs  considerably  exceeds  the  error  estimates.  The 
reason  for  this  is  unknown. 

The  theories  chosen  for  comparison  are  those  of  Stewartson-Wedemeyer 1 , 
Kitchens,  et  al.9,  and  Gerber7.  The  first  two  yield  only  eigenvalues  for  the 
comparison.  A  basic  assumption  of  the  Stewartson-Wedemeyer  theory  is  that 
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near  resonance  of  an  eigenmode,  the  response  of  this  mode  dominates  the 
response  of  the  liquid.  This  implies  that  if  the  liquid  is  forced  by  steady 
coning,  the  peak  of  Cp  should  occur  when  w  is  equal  to  an  eigenfrequency. 

Thus,  the  proper  comparison  is  between  the  computed  eigenfrequency  and  the 
value  of  w  at  which  the  peak  of  Cp  occurs.  The  Gerber  theory  yields  Cp  as  a 

function  of  w,  and  should  properly  be  compared  at  all  frequencies. 

The  comparisons  are  presented  in  Figure  10(a)  to  10(g).  In  each  case, 
the  vertical  dashed  line  indicates  the  eigenfrequency;  on  the  scale  of  the 
plots,  the  Stewartson-Wedemeyer  and  Kitchens,  et  al.,  predictions  do  not 
differ  visibly.  The  two  vertical  dotted  lines  show  the  eigenfrequencies 
computed  for  cylinders  with  c/a  incremented  and  decremented  by  its  estimated 
error.  For  the  partially-filled  cases,  the  dotted  lines  also  incorporate  the 
estimated  error  in  fill  ratio.  The  solid  curve  in  each  figure  is  the  Gerber 
prediction;  no  curve  is  shown  for  the  partially-filled  cases,  since  the  theory 
does  not  apply  there. 

The  dependence  of  the  predicted  and  observed  frequencies  of  peak  response 
on  Reynolds  number  is  shown  in  Figures  11(a)  and  (b).  In  these  figures  the 
size  of  the  symbols  roughly  corresponds  to  the  effect  of  aspect  ratio  error. 
In  each  ease,  Stewartson-Wedemeyer  prediction  seems  to  fall  away  from  the  data 
at  lower  Reynolds  number;  for  Re  =  5  *  103,  c/a  =  3.1481,  the  error  is  about 
eight  percent.  The  Gerber  prediction  matches  the  data  much  better,  the  worst 
error  being  about  two  percent. 

Figure  12  shows  the  value  of  Cp  at  its  peak,  plotted  against  Re.  Again, 

the  size  of  the  symbols  indicates  the  error.  The  solid  lines  are  the  Gerber 
prediction.  The  data  indicate  a  less  rapid  decrease  of  peak  Cp  with  decreas¬ 
ing  Re  than  does  the  theory. 

All  of  the  data  presented  so  far  were  taken  using  a  transducer  at  r/a  = 
0.668.  A  second  transducer  was  also  mounted  in  the  end  wall,  at  r/a  =  0.436. 
For  various  reasons,  data  were  taken  using  both  transducers  only  about  one- 
third  of  the  time.  The  data  taken  with  the  second  transducer  showed  the  same 
behavior  as  those  taken  at  r/a  =  .668,  as  to  linear  and  nonlinear  behavior  and 
location  of  peak  response.  However,  Cp  at  r/a  =  .436  was  consistently  less 

than  at  .668,  which  is  an  indication  of  the  form  of  the  radial  structure  of 
the  fluctuating  pressure  field.  A  measure  of  the  variation  between  radial 
locations  is  the  ratio  Cp  (r/a  =  .436) /Cp  (r/a  =  .668).  This  ratio  is  plotted 

against  w  in  Figures  13(a)  to  (e)  for  those  cases  in  which  data  are 
available.  Note  that  in  these  figures  all  available  data  are  plotted, 
including  those  in  the  nonlinear  regime.  The  vertical  dashed  line  shows  the 
observed  frequency  of  peak  Cp  for  the  linear  regime.  The  solid  line  indicates 

the  prediction  of  the  Gerber  theory. 


D.  The  Nonlinear  Regime 

The  location  and  nature  of  this  regime  depended  on  Re,  c/a,  and  fill 
ratio.  Figures  14(a)  to  (g)  show  Cp  plotted  against  w  for  all  available  data, 

at  r/a  =  .668.  The  different  symbols  indicate  different  coning  angles. 

Away  from  the  peak  of  Cp,  most  of  the  data  in  each  plot  seem  to  lie  along 

the  same  curve,  showing  linear  behavior.  For  sufficiently  large  e,  the  curve 
seems  to  be  chopped  off  near  the  location  of  the  linear  peak.  This  effect 
would  not  be  well-modeled  by  a  simple  shift  in  location  of  the  peak  response, 
so  the  theory  of  Scott10  is  not  supported.  For  still  larger  e,  the  behavior 
becomes  nonstationary,  indicated  as  before  by  vertical  lines  connecting  data 
points. 

A  measure  of  nonlinearity  is  the  change  in  Cp  at  the  frequency  at  which 

Cp  peaks  in  the  linear  regime  —  call  this  value  Cp*.  There  is  some 

arbitrariness  in  this  measure,  since  the  nonlinear  regime  Cp  sometimes  peaks 

at  a  different  frequency,  but  it  serves  to  show  the  onset  of  different  types 
of  behavior.  Figures  15(a)  to  (g)  show  Cp*  plotted  against  e.  The  error  bars 

now  indicate  uncertainty  in  estimating  Cp*.  and  are  based  on  the  observed 

scatter  of  the  data  rather  than  the  estimated  uncertainty  of  the  measurement 
of  Cp.  The  symbols  connected  by  vertical  arrows  indicate  nonstationary 

behavior. 

The  onsets  of  nonlinear  and  nonstationary  behavior  can  be  seen  to  be 
strongly  influenced  by  Re;  for  Re  -  5  *  103,  nonstationary  behavior  was  not 
observed  even  at  the  maximum  coning  angle  attained,  e  =  0.035.  Figure  16 
shows  the  values  of  e  at  which  transitions  to  nonlinear  behavior  occurred,  as 
functions  of  Re.  Figure  17  shows  a  similar  plot  of  the  transitions  to  non¬ 
stationary  behavior.  Note  that  in  Figure  17,  only  a  lower  bound  is  indicated 
for  c/a  =  3.1481,  Re  =  5  x  io3,  since  nonstationary  behavior  was  not  observed 
for  this  case. 


IV.  DISCUSSION  AND  CONCLUSIONS 
A.  Validity  of  Linear  Theory 

The  Stewartson-Wedemeyer 5  and  Kitchens,  et  al.9  eigenfrequencies  approach 
the  Gerber7  frequency  of  peak  response  as  Re  becomes  large  (Figures  11(a)  and 
(b)).  The  data  agree  with  all  of  these  theoretical  results,  to  within  1.5%  in 
the  worst  case. 

As  Re  becomes  small,  however,  the  location  of  the  peak  of  the  Gerber 
theoretical  curve  diverges  from  the  eigenfrequency.  The  data  agree  better 
with  the  Gerber  theory  (within  1.5%)  than  with  the  eigenfrequencies  (greatest 
difference  of  9%).  This  divergence  of  the  frequency  of  peak  response  from  the 
location  of  the  S-W  eigenfrequency  probably  reflects  the  increasing  response 
of  the  other  inertial  modes.  Since  the  Gerber  theory  computes  the  response 
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of  these  other  modes,  better  agreement  occurs.  At  low  Re,  the  response  of  the 
liquid  may  be  controlled  by  several  modes. 


The  predictions  of  the  Gerber  theory  of  the  maximum  response  (Figure  12) 
differ  from  the  data  by,  at  most,  131.  The  maximum  difference  occurs  at  the 
largest  value  of  Re,  which  is  not  consistent  with  an  error  due  to  improperly 
modeled  viscous  effects.  The  trend  with  decreasing  Re  also  differs  from  the 

data,  with  the  predicted  peak  Cp  varying  approximately  with  Re1/2,  and  the 

data  varying  approximately  with  Re**3  for  c/a  =  3.1481,  and  with  Re*3^  for  c/a 
=  1.0509. 

The  limited  data  available  indicate  some  difference  between  the  predicted 
and  observed  radial  structure  of  the  two  modes  observed.  There  is  no 

indication  that  this  is  the  cause  of  the  observed  differences  in  peak  Cp, 

however,  since  the  discrepancies  in  radial  structure  at  the  frequency  of  peak 
Cp  have  no  apparent  relation  to  the  discrepancies  of  Cp  itself.  The 

differences  do  serve  to  qualify  the  remark  made  in  Section  I IA  of  this  report 
regarding  the  location  of  the  pressure  transducers.  If  the  observed  spatial 
structure  differs  from  the  predicted  structure,  then  it  is  conceivable  that 
the  extent  of  agreement  found  here  is  not  representative  of  the  whole  flow. 
Only  measurements  at  a  large  number  of  points  would  resolve  this  issue. 

B.  Applicability  of  Linear  Theory 

The  onset  of  noticeably  nonlinear  behavior  occurs  at  larger  c  for  smaller 
Re,  but  even  at  Re  =  5  x  103,  the  smallest  value  of  Re  used  in  this  experi¬ 
ment,  the  transition  occurred  at  e  =  4  x  10" 3  rad.,  which  is  much  smaller  than 
precession  angles  typical  of  artillery  projectiles.  Since  the  existing 

theories  are  not  applicable  for  small  Re  (less  than,  say,  103) ,  their  direct 

applicability  to  any  artillery  projectiles  is  questionable.  This  does  not 
mean  that  they  are  useless,  however,  since  they  may  serve  as  the  basis  for 
theories  (perhaps  semi -empirical)  which  can  account  for  nonlinear  effects. 

C.  Nonlinear  Behavior 


Based  on  work  cited  in  the  Introduction  17'18,  it  is  reasonable  to 
hypothesize  that  the  onset  of  nonstationary  behavior  is  associated  with  an 
instability  of  a  secondary  flow  induced  by  the  forced  inertial  mode.  If  this 
is  so,  then  there  is  some  time  scale  associated  with  growth  of  the  mode  to  a 
point  of  instability.  Since  free  gyrostat  experiments  are  inherently  time- 
dependent,  this  may  account  for  the  fact  that  in  free  gyrostat  experiments 
departure  from  linear  behavior  was  observed  at  larger  e  than  in  the  present 
experiment. 


V.  RECOMMENDATIONS  FOR  FUTURE  WORK 

The  most  obvious  question  raised  by  this  experiment  concerns  the  predict¬ 
ability  of  the  nonlinear  and  nonstationary  behavior.  If  the  time  scale  of  the 
fluid  instability  mentioned  in  the  previous  section  is  longer  than  the  flight 
time  of  a  typical  projectile,  then  the  proper  theoretical  approach  is  to  model 
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the  non-steady  nonlinear  evolution  of  the  flow  before  It  becomes  unstable.  If 
the  time  scale  is  shorter  than  the  flight  time,  then  the  details  of  the 
approach  to  instability  are  irrelevant,  and  instead  some  attempt  must  be  made 
to  model  the  flow  after  the  Instability  occurs,  and  the  flow  becomes 
(probably)  turbulent.  Thus,  the  first  experimental  task  is  to  determine  the 
time  scale. 

If  it  turns  out  that  the  time  scale  is  long,  and  that  the  evolution  of 
the  flow  before  it  becomes  unstable  is  important  to  the  question  of  projectile 
stability,  then  two  Issues  mist  be  addressed:  first,  whether  the  spatial 
variation  of  the  flow  field  is  essentially  as  described  by  linear  theory,  with 
perhaps  small  departures,  and  second,  whether  the  Kitchens,  et  al.,  theory 
correctly  describes  the  effect  of  non-rigid  rotation  on  the  resonant 
frequencies.  If,  on  the  other  hand,  the  post- instability  flow  is  of  greatest 
importance,  then  it  must  be  determined  whether  the  flow  is  amenable  to  some 
statistical  or  other  description  which  allow  the  moment  exerted  on  the 
cylinder  to  be  specified  sufficiently  well  to  predict  system  response. 

The  present  experiment,  and  the  ones  suggested  for  the  future,  are  felt 
to  be  relevant  to  projectiles  and  gyrostats,  despite  the  fact  that  e  is  held 
constant.  However,  this  should  be  demonstrated  directly  by  measuring  unsteady 
pressure  in  a  free  gyrostat  or  in  projectiles. 
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CYLINDER  ROTATION  AXIS 


Figure  2.  Glmbal  Case 
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Figure  3.  Adjustable  Cam. 
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Figure  5.  Telemetry  Link. 
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Figure  7.  Overlaid  Spectra  Taken  at  Different  Coning 
Frequencies. 


Figure  9.  Variation  of  Cp  with  Coning  Frequency, 

c/a  =  3.1481,  f  =  0.92,  Re  =  4  x  105 
A  -  e  =  1.85  x  10"H;  □  -  e  ■  3.51  x  10" *♦ 
x  -  E  =  7.78  x  10-4;  0  -  e  -  1.97  x  10" 3 
O  -  e  =  3.91  x  10- 3. 
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Figure  10.  Variation  of  Pressure  Coefficient  with  Coning  Frequency 
in  the  Linear  Regime.  Transducer  at  r/a  =  0.668. 

Dashed  vertical  line  is  eigenfrequency  at  nominal  c/a. 
Ref.  5,  6,  and  9;  dotted  vertical  lines  show  the  effect 
on  eigenfrequency  of  experimental  error  in  c/a.  Solid 
curve  is  prediction  of  Ref.  7. 

a.  c/a  =  3.1481,  f  =  1.00,  Re  =  5  x  105. 
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Figure  lOg.  c/a  =  3.1481,  f  =  0.92,  Re  =  8 


Figure  11.  Variation  of  Frequency  of  Peak  Response  with 
Reynolds  Number.  Transducer  at  r/a  =  0.668. 
Solid  Lines  are  Stewartson-Wedemeyer 
eigenfrequency;  dashed  lines  are  Gerber 
frequency  for  peak  response. 


Figure  12.  Variation  of  Amplitude  of  Peak  Response  with 
Reynolds  Number.  Transducer  at  0.668. 

Solid  lines  are  Gerber  peak  response, 
f  =  1.00;  0  -  c/a  =  3.1481;  a  -  c/a  «  1.050 
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Figure  13.  Ratio  of  Amplitude  at  r/a  =  0.436  to  Amplitude  at 
r/a  =  0.668.  Solid  line  is  Gerber  theory. 

a.  c/a  =  3.1481,  f  «  1.00,  Re  *  5  x  105 

A  -  e  =  1.66  x  10"4;  □  -  e  =  3.79  x  10~4; 

X  -  e  =  5.93  x  10'4;  0  -  e  =  8.17  x  10"4; 

V  -  e  =  3.84  x  10-3. 
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Figure  13b.  c/a  =  3.1481,  f  =  1.00,  Re  =  1  *  105 

A  -  e  =  1 . 75  *  10'4;  □  -  e  =  3.79  x  10'4. 
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Figure  14.  Variation  of  Pressure  Coefficient  with  Coning  Frequency 
in  the  Linear/Honlinear/Nonstationary  Regimes  Transducer 
at  r/a  =  0.668.  Symbols  connected  by  vertical  lines 


represent  nonstationary  behavior. 

a.  c/a  =  3.1481,  f  =  1.00,  Re  =  5  *  10s 

A  -  E  =  2  *  10’4;  □  -  e  =  3.79  x  10'4; 

X  -  E  =  6  x  10"4;  </  -  e  =  8.17  x  10'4; 

O  -  e  =  1.85  x  10"3;  V  -  e  =  3.84  x  10"3. 


Figure  14b.  c/a  =  3.1481,  f  =  1.00,  Re  =  1  *  105 

A  -  e  =  1.75  x  10"4;  □  -  e  =  4  x  10'4 

X  -  e  =  9.23  x  10"4;  0  -  e  =  3.68  x  1 


Figure  15.  Dependence  of  Peak  Amplitude  on  Coning  Angle, Symbols 
Connected  by  Vertical  Arrows  Represent  Nonstationary 
Behavior.  Error  bars  indicate  uncertainty. 


a.  c/a  =  3.1481,  f  =  1.00,  Re  =  5  *  105. 


Figure  159-  c/a  =  3.1481,  f  =  0.92,  Re  =  8  *  10\ 
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Figure  16.  Coning  Angle  at  Which  Transition  to  Nonlinear  Behavior  Occurred 
O  -  c/a  =  1.0509,  f  =  1.00;  a  -  c/a  =  3.1481, 
f  =  1.00;  □  -  c/a  =  3.1481,  f  =  0.92. 


Figure  17.  Coning  Angle  at  Which  Transition  to  Nonstationary 
Behavior  Occurred 


LIST  OF  SYMBOLS 
(Continued) 


coning  frequency 


70 


APPENDIX  A:  PROPERTIES  OF  LIQUIDS 


The  liquids  used  in  this  experiment  were  Dow  Corning  200  Silicon  Fluid, 
of  lcS,  5cS,  and  lOOcS  nominal  kinematic  viscosities.  Relevant  properties, 
extracted  from  the  manufacturer's  literature  are  shown  in  Table  A-l. 


Table  A-l.  Properties  of  Liquids 


Viscosity  at  25°  C 
(nominal ) 

Specific  Gravity 
at  25°  C 

Coefficient  of 
Thermal  Expansion  (B) 
(cc/cc/°C) 

lcS 

0.818 

.00134 

5cS 

0.940 

.00105 

lOOcS 

0.960 

.00096 

Also  in  the  manufacturer's  literature  are  charts  showing  the  variation  of 
viscosity  with  temperature,  for  the  5cS  and  lOOcS  fluids.  For  convenience  in 
computer  data  reduction,  a  simple  function  of  the  form  v(T)  =  exp  (A+BT)  was 
fit  to  several  points  read  from  each  of  the  two  manufacturer's  curves,  and  was 
used  to  determine  the  viscosity  of  each  liquid  as  a  function  of  room  temper¬ 
ature  at  the  time  of  each  run.  For  the  range  of  temperature  encountered  in 
the  experiment  (19.6°  C  -  27.8°  C)  such  a  function  fit  the  points  read  from 
the  charts  to  within  one  percent  greatest  error,  which  probably  is  less  than 
the  error  in  reading  the  chart. 

The  U.S.  Army  ARRADCOM  Chemical  Systems  Laboratory  measured  the  viscosity 
of  a  sample  of  the  lcS  fluid  used  in  the  experiment,  finding  it  to  be  1.044  cS 
at  25°  C  and  0.932  cS  at  35°  C.  The  same  simple  exponential  used  above  was 
fit  to  these  points,  and  was  used  in  the  same  fashion.  The  Laboratory  also 
measured  the  density  at  25°  C,  obtaining  0.8149  g/ml  as  the  result.  The 
coefficients  used  in  the  viscosity  curves  are  shown  in  Table  A-2.  To  estimate 
the  variation  of  density  with  temperature,  it  was  assumed  that  the  coefficient 
of  thermal  expansion  B  remained  constant  with  temperature.  With  this  assump¬ 
tion  it  is  easy  to  show  that 

p(T)  =  p(25°  C)  •  [1  -  p(25°  C)*3-  (T-25°  C)] 
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Table  A-2.  Coefficients  in  Expression  v(T)  =  exp  (A+BT) 


Fluid 

(Nominal  Viscosity) 

A 

B 

lcS 

.35402 

-.012339 

5cS 

2.1832 

-.02122 

lOOcS 

5.1312 

-.020655 
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APPENDIX  B:  CALIBRATIONS 


1.  Pressure  Transducers.  The  pressure  transducers  were  calibrated  in  two 
stages: 


a.  Three  transducers,  called  "callbratio.;"  transducers,  were  calibrated 
against  a  Consolidated  Electrodynamics  Type  6-2C1  Primary  Pressure  Standard. 
This  calibration  was  performed  twice,  once  before  and  once  after  the  comple¬ 
tion  of  the  experiments.  The  resulting  pressure  versus  voltage  data  were  fit 
with  three  straight  lines,  one  for  each  transducer.  The  slope  of  the  lines 
changed  about  one  percent  (worst  case)  between  the  two  calibrations. 

b.  The  two  transducers  in  the  cylinder,  called  "test"  transducers,  were 
calibrated  against  the  other  three  transducers.  This  calibration  was  done 
four  times  during  the  experiments.  The  slopes  of  the  resulting  calibration 
curves  changed  by  no  more  than  0.4%  (worst  case)  between  calibrations  of  the 
same  transducer. 

Including  uncertainties  in  both  calibrations,  the  final  calibration 
constant  Kp  was  probably  accurate  to  about  one  percent. 

2.  Tunable  Discriminators.  The  EMR  model  4142  tunable  discriminators  had 
adjustments  for  deviation,  coarse  center  frequency,  offset  (fine  center  fre¬ 
quency),  type  of  low-pass  filter,  filter  cutoff  frequency,  and  bandedge  volt¬ 
age.  The  coarse  center  frequency  was  always  set  to  the  VCO  nominal  center 
frequency  (40  kHz  or  52.5  KHz).  The  low  pass  filter  used  was  of  constant- 
delay  type,  with  a  cutoff  frequency  of  1  KHz.  The  deviation  was  set  to  7.5%, 
the  same  as  the  nominal  VCO  deviation. 

The  input  of  each  VCO  could  be  set  to  the  transducer  amplifier  output 
(the  data  signal),  to  a  D.C.  reference  signal  Vp,  which  was  approximately  2.5 

volts  above  VCO  reference,  or  to  VCO  reference,  VR.  The  discriminator  offset 

was  adjusted  so  that  when  the  VCO  input  was  set  to  Vq,  the  discriminator 

output  was  (V0  -  2.50)  +  .005  volts.  The  discriminator  bandedge  voltage  was 

adjusted  so  that  when  the  VCO  input  was  set  to  VR,  the  discriminator  output 

was  -2.50  +  .005  volts.  After  these  adjustments  were  performed,  the  discrim¬ 
inator  output  was  equal  to  the  VCO  input  minus  2.50  +  .03  volts,  at  any  VCO 
input  voltage  from  zero  to  five  volts. 

3.  System  Transfer  Function.  The  calibration  of  transducer  pressure  to 
discriminator  output  voltage  was  conceptually  divided  into  two  parts:  the 
conversion  of  pressure  to  voltage  by  the  transducer,  represented  by  the  trans¬ 
ducer  calibration  coefficient,  Kp,  and  the  various  voltage  level  changes, 

voltage-to-frequency  conversions  and  their  inverses  accomplished  by  the  trans¬ 
ducer  amplifiers  and  the  telemetry  link,  represented  by  the  system  transfer 
function,  H(A).  The  measurement  of  the  transfer  function  was  accomplished 
with  the  Hewlett-Packard  3582A  spectrum  analyzer  (which  was  also  used  to 
collect  the  experimental  data).  To  perform  the  measurement,  a  pseudo-random 
noise  source,  available  as  an  output  from  the  analyzer,  was  connected  to  a 
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transducer  amplifier  input,  in  parallel  with  the  transducer,  and  also  to  one 
input  channel  of  the  analyzer  (Figure  B-l).  The  output  of  the  corresponding 
discriminator  was  connected  to  the  other  analyzer  input,  and  the  analyzer  was 
adjusted  so  as  to  perform  exponential  averaging  over  16  sample  spectra,  then 
compute  the  transfer  function  from  one  input  signal  to  the  other  using  the 
averaged  spectra  (see  the  manufacturer's  literature22  for  details).  This 
process  was  repeated  for  the  other  channel,  and  the  results  were  stored  in  the 
HP  9845B  computer  which  was  connected  to  the  analyzer.  Note  that  the  computed 
transfer  function  was  complex;  both  magnitude  and  phase  were  stored. 

4.  Determination  of  Fluctuating  Pressure  Amplitude.  The  discriminator  output 
was  a  mixture  of  signals  of  the  form  Vd(x,t)  =  V(x)  cos  (xt  +  6), 
corresponding  to  fluctuating  pressure  signals  p(X,t)  =  p(x)  cos  (Xt  +  $).  The 
spectrum  analyzer  presented  V(X);  the  amplitude  of  a  particular  spectral 
component  of  the  fluctuating  pressure  was  then 

PU)  «  Kp  V(X)  /  | H ( X ) | , 

where  Kp  =  dp/dV,  the  slope  of  the  pressure  transducer  calibration  and 
| H (x ) |  is  the  magnitude  of  the  telemetry  system  transfer  function. 
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